RECURSIVE FORMULAS FOR SUMS OF SQUARES AND SUMS 
OF TRIANGULAR NUMBERS 
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Abstract. We prove recursive formulas for sums of squares and sums of tri- 
angular numbers in terms of sums of divisors functions and we give a variety of 
consequences of these formulas. Intermediate applications include statements 
about positivity of the coefficients of some infinite products. 



1. Introduction 

The purpose of this paper is to give inductive formulas for sums of squares, sums 
of triangular numbers, and mixed sums of squares and triangular numbers. Our 
main tool is the following theorem. Let N = {1, 2,3,.. .} and let No = {0, 1,2,.. .}. 

Theorem 1. 3J Let a e N. Let A 1: A 2 , . . . A a C N and A = {A u A 2 , . . . , A a ) 

and let fa : Ai — > C for i = 1,2,..., a be arithmetic functions and let f — 
(/i,/ 2 ,...,/«). U both 

a oo 

Fa{x) = J] II 0- - x")-^ = 5>A,/(n)a;« 

i=l ngA, n=0 



EE 

i=l n£A, 



converge absolutely and represent analytic functions in the unit disk \x\ < 1, then 

n / a. \ 

npAj(n) = ( PA,f(n - k) ^ fuA t i k ) 
where paj{0) = 1 and 



/ j i v '"J • 

k=l \ i=l 



AaM= E Md)- 

d\k 
deA z 

We note that Theorem [T] for the special case a = 1, has been given in Apostol 
[I] and in Robbins in [6] to give formulas relating arithmetic functions to sums of 
divisors functions. The authors' key argument is that generating functions for the 
arithmetic functions they considered have the form of infinite products ranging over 
a single set of natural numbers. In our previous work [3] we used Theorem [1] to deal 
with arithmetic functions whose generating functions involve finitely many infinite 
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products ranging over different sets of natural numbers. We derived a variety of 
inductive formulas for such functions. In the present note we shall continue employ- 
ing Theorem [T] to deduce inductive formulas for sums of squares, sums of triangular 
numbers, and mixed sums of squares and triangular numbers, see Section 2. We 
notice at this point that there is a large literature on sums of squares and sums 
of triangular numbers. The interested reader is referred to the classical volumes of 
Dickson [2] and the recent book of Williams [7] along with their references. With 
an essential help of identities given in Williams in [7], we will give a variety of 
consequences of our recursive formulas, see Section 3. Further, we shall prove that 
the coefficients of some infinite products are all positive, see Section 4. By way of 
example, we will prove the following results: 

(1) 

^2(a(n)-Aa(n/4))xA (j^ (<r*(n) - 4o>/2)) xA = 

OO 

E (n (o-(rc) - 4a(n/4)) - (a*(n) ~ Aa*(n/2)))x n , 

n=l 

see Corollary [1] below. Here <r(n) and <J*(n) are as in Definitions [1] and [3] below. 

(2) If p = 1 mod 4 is prime, then 

E (**(p - 3) 4 ^(^)) = p - L 

3=1 ^ ' 

refer to Theorem [6] below. Here ri (n) is the number of representations of n as a 
sum of two squares. 

(3) If p and 4p + 1 are primes, then 

X>C?) Up-^-M^-)) =-i, 

see Theorem [8] below. Here t2{n) is the number of representations of n as a sum of 
two triangular numbers. 

(4) The following is a direct consequence of Theorem [11] below. If x G C with 
\x\ < 1 and a £ N, then the coefficients of the infinite series 

oo oo 

Y[(l -x n )- a {l -x 3n ) a , Y[(l-x n )- a {l-x 3n - 1 ) a , 

n—1 n— 1 

and 

oo 

Y[(l -x n )- a (l -x 3n ) a (l -x 3n ~ 1 ) a 
n=l 



are all positive. 

Our identities involve sums of divisors functions which are introduced in the 
following three definitions. 



Definition 1. Let the function a be defined on Q as follows: <j(0) = 1, if q £ ^ \no, 
then a(q) = , and if q £ N, then 

a(q) = J2d. 

d\q 



RECURSIVE FORMULAS FOR SUMS OF SQUARES.. 



Definition 2. Let n,r £ No, let m £ N, and let 



d\n 
d=r mod m 

If m = 2 and r = 1 we often write cr°(n) rather than o-i^{n) and if m = 2 and 
7' = we often write a E (n) rather than 170,2 (n)- 

We have the following basic facts: for m,n £ N 

(1) <r(n) — <r°(n) + o~ E {n), cr-m,2m{ n ) — mo~°(n/m), and cro. m (n) = mo~(n/m). 

Definition 3. Let the function a* be defined on Q as follows: if q £ Q \ N, then 
<r*(q) = and if q £ N, then 

a*(q)= d - 

d\q, ^ odd 

By [71 Theorem 3.4] we have for all n £ N that 

(2) a*(n)=a(n)-a(n/2). 

Further we will need the following two identities due to Jacobi and Gauss respec- 
tively: 

00 00 

(3) Wil-x^l + x 2 "- 1 ) 2 = x " 2 - 

n—1 7i— — 00 

(4) na-a^xi-^r^f;^. 

n— 1 n— 

Identities © and (0} can be found for instance in Hardy and Wright [I]. 

2. General formulas 

In this section we give inductive formulas for sums of squares, sums of triangular 
numbers, and mixed sums of squares and triangular numbers. 

Definition 4. Let k £ No and let the function r k (n) be defined on No as follows: 

r k (n) = #{(xi,x 2 , . . . ,x k ) £Z k : n = x\ + x\ + . . . + x 2 k }. 

Theorem 2. If n £ N, then 

nr k (n)=2k I a*{n) - 4a*(n/2) + ]T r k (j) ( a * (n - j) - 4a* (^-) 



3=1 x 

Proof. Taking k powers in identity (J3)) we have 

/ 00 \»= 

JJ (1-^(1 + ^-1^= J2 ^ =l + ^r fc (n)x", 

n—1 \n— — 00 / n—1 

or equivalently, 

00 00 

JJ(1 - x 2n ) fc (l - x in - 2 ) 2k {l - x 2n - 1 )~ 2k = l + Yl r k {n)x n . 
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Let in Theorem [U A x = 2N and fi(n) = -kn, A 2 = 4N - 2 and f 2 (n) = -2kn, 
and A 3 = 2N - 1 and f 3 (n) = 2kn. Then for n > 1 

nr^n) = -ka E (n) - 2ka 2 ,A(n) + 2ka°(n)+ 

(5) n_1 

^2r k (j)(-ka E (n-j)-2ka 2 An-j)+2ka (n-j)). 

j'=i 

Identities ([T]) and formula ([2]) yield 

-a E (n) - 2a 2A (n) + 2a°(n) = 2 (<r(n) - <r(n/2)) - 8 (cr(n/2) - a(n/4)) 

= 2<r*(n) - 8<T*(n/2). 

Now put in formula ([5]) to conclude the desired identity. □ 
Definition 5. Let k 6 No and let the function t k be defined on No as follows: 

t*(n) = # | (xx, x 2 , ■ ■ ■ , x k ) € N k : n = £ __v_|±_) | . 

Theorem 3. If n € No, i/ien 



ra-1 / _ o \ 

ntfc(n) = *5D**0') (^(n-jJ-M-o-^)) • 
Proof. Taking k powers in identity we have 

oo / oo \ ^ oo 

n (i - * 2n ) k (i - * 2n - l T k - E ^ = E **(»)^. 

n=l \n=0 / n=0 

and the result follows by TheoremQ] applied to A\ — 2N, = —kn, A 2 = 2N— 1, 

and f 2 (n) = kn. □ 

Definition 6. Let k, I € N and let the function tifc^ be defined on No as follows: 

uu(n) = #Ux 1 ,...,x k ,y 1 ,...,y l ) e Z k x N : n = E^+E ^ ± ^ V 
t i=i i=i J 

Theorem 4. J/ k,l,n G N, then 

nu kt i(n) = {2k + l)a(n) - 2(5fc + 2/)<r(n/2) + 8fca(n/4)+ 



E «fc,l(7) ( 2fc + - 3) - 2(5fc + 2l)a{^l) + 8fca(^ 

j'=i V 

Proof. Clearly 

oo oo 

JJ(1 - X 2n ) k+l (l + a; 2 "- 1 ) 2 ^! - x 2 "- 1 )- 1 = 1 + E «W 

n— 1 n—1 

that is, 

oo oo 

JJ(1- z 2n ) fe+ '(l - x 4 "- 2 ) 2fe (l - x 2 "^ 1 )- 2 ^! - x 2 ^ 1 )- 1 = 1 + E «fc,i(")3 
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equivalently, 

oo oo 

JJ(1 - x 2n ) 3k+l (l - x 4n )- 2k {\ - a**- 1 )-**- 1 = l + Y, u k ,i(n)x n , 

n—l n—1 

which by Theorem [Tj gives the desired identity. □ 

3. Sums of two, four, and eight squares 

The following formulas for ^(n), r^n), and fg(n) are well-known and can be 
found for instance in Williams [7]. If n G N, then 

d\n ^ 

(6) r 4 (n) = 8<r(n) - 32a(n/4), 

r 8 (n) = 16(-l)"^(-l) d d 3 

d\n 

where 

1, if d = 1 mod 4, 
— 1, if d = 3 mod 4, 
0, otherwise. 
We start by evaluating a convolution. 

Theorem 5. If n G N, t/ien 

E( CT 0')-Mf)) (^(n-i)-4 C r*(^, 

n(a(n)-4a(^))-(a*(n)-4a*(^)). 
Proof. By Theorem [2] 

n— 1 / . v 

(7) nr 4 (n)=8 ( 7*(n)-32 ( 7*(n/2) + Er4(i) I 8a*(n - j) - 32a*(^) J , 

which combined with (|6]) for r4(n) yields the result. □ 
Corollary 1. 

(X>(n)-4a(n/4))^ (j^(a*(n) - 4a*(n/2)) xA = 

OO 

(ncr(n) - 4na(n/4) - a*{n) + 4a*(n/2))x n . 

n=l 

Proof. Immediate from Theorem [5] □ 

Theorem 6. Let p be prime, 
(a) If p = 1 mod 4, i/jen 

£ ^0') (<^*(p - i) - ^(^)) = p - 1. 
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(b) If p = 3 mod 4, then 
p-i 



E ^00 ( °*{P - J) 4a*(^) ) = -p - 1. 



2 

Proof, by Theorem [5] we have 



(8) nr 2 (n) = 4a* (n) - 16a* (n/2) + E r a (?) (^(n - j) - 16a*(^)) . 

Further, if p is prime, then clearly a*(p) = 1 + p and a*(p/2) = 0. Moreover, by 
virtue of formulas ([6]) we have that ^(p) = for p = 3 mod 4 and r2(p) = 8 for 
p = 1 mod 4. Now combine these facts with identity ([8]) to conclude parts (a) and 
(b). □ 

Corollary 2. Let (p,p + 2) be a twin-prime, 
(a) If p = 1 mod 4 ; i/ien 
P+i 



E ^(J) (**(P + 2 - i) - 4a*( P + 2 J )) 
j= i V / 

p- 1 / 

tP-3 



-4-E^(j) U*(p-j)-4a*(^ 
f&J If p = 3 mod 4, i/ien 

E r * (?) ^> + 2 - i) - 4a*( P + 2 J 

Proof. Immediate from Theorem [B] □ 
Theorem 7. Lei p &e an odd prime. Then 

(a) i>(j) (Vb-j)-4a*(^)) = p 2 -l. 



p-i 



(6) £r 8 (j) a*(p-j)-4a*(^) = P 4 -1- 

Proof, (a) As mentioned before a*(p) = l+p and a*(p/2) = 0. By identity (|6|) we 
have r4(p) = 8(1 + p). Putting in formula ([7]) we get 

8p(l + p) = 8(p + 1) + E »-4(3) (8a> - i) - 32a*(^) > ) , 

which gives the desired identity, 
(b) By Theorem [3 

(9) nr 8 (n) = 16 ( a*(n) - 4a*(n/2) + £ r S (j) (a*(n - j) - 4a*(^)) 
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If p is an odd prime, then by formulas ([6]) we have rg(p) = 16(1 + p 3 ), which 

combined with identity ([9]) gives the result. □ 

Corollary 3. Let p be an odd prime. Then 



(a) h + ^r 2 (j)(<7*(p-i)-4a*(^) 
1 + X> 4 (j) [a*{p-j)-Aa*{ P -^-) 

(b) |l + ^ 4 (j)^(p^VV(^) 

1+E r *« (o-*(p-j)-4a*(^±) 

Proof. Immediate from Theorems |5] and [7J □ 

4. Sums of two, four, and six triangular numbers 

The following formula for t2(n), ti(n), and te(n) can be found in Ono, Robins, 
and Wahl [5] and in Williams [7]. 

d|4n+l V 

(10) t 4 (n) = <r(2n + 1), 

-4 



d|4n+3 

Theorem 8. If p and 4p+ 1 are prime numbers, then 

X>(j) (<7(p-i)-M^)) = -i. 

Proof. Suppose that p and 4p+ 1 arc primes. Clearly <j(p) = 1 +p, and by formulas 
(ITOl) we have i 2 (p) = 2. Further by Theorem [3] 

p-1 / 

(11) pt 2 (p) = 2a(p) - 8<r(p/2) + E *2(j) MP ~ j) ~ M^ 2 ) 

i=i V 

Thus 

2p = 2(1 +p) + 2^t 2 (j) ff (p -j) - 4 CT (^— i) 



Theorem 9. If 2n + 1 is prime, then 



□ 



EMi)(^-i)-M^))=^±ll (=Ej). 



j=0 ^ ' i=i 



8 MOHAMED EL BACHRAOUI 

Proof. Assume that 2n + l is prime. Then by identities (fT0|) we have ti{n) = 2n + 2. 
Moreover, from Theorem [3] we get 

n—l / 

(12) nU(n) = 4 £ t 4 (j) ( a(n - j) - 

Combining these two identities proves the result. □ 
Theorem 10. If 4n + 3 is prime, then 

E" . ,.J , . \ , ,n — jA n(n + l)(2n + l) . . 2 . 
*eO)K»-j)-M-2^)J =— 7T " (=2^)- 

Proof. If 4n + 3 is prime, then by (|10[) we have 

*e(«) = ~(1 - (4n + 3) 2 ) = (n + l)(2n + 1). 
o 



Next by Theorem [3] we have 



n-l 

n — j 



nt G {n) = 6 ^ t 6 (j) ( cr(n - j) - 4cr( 



2 

Using the previous two formulas we obtain the result. □ 

5. Facts on some infinite products 

Throughout this section we suppose that x is a complex number such that \x\ < 1. 

Theorem 11. Let a, b, n € N and let I be a nonempty subset of {0,1,2, ... ,b — 2} . 
Then the coefficients of the infinite product 

oo 

n=i iei 

are all positive. 
Proof. Write 

oo oo 

n lit 1 - xn y a ( i - x hn ~T = A (n)x n . 

n—l iei n—0 

We show by induction that A{n) = for all n £ N . Clearly A(0) = 1 > 0. Now 
assume that the assertion is true for j = 0, 1, . . . , n — 1. By Theorem [T] we have 

n-l / \ 

nA(n) = aa(n) - a^Vj^n) + ^ A(j) I aa(n - j) - a^cr ijb (n - j) J . 

iei 3=1 V iei / 

But 

6-1 

aer(n) - a^er j; . b (n) = a} a ijb (n) - a^Vj :b (n) > ai :b (n) > 0. 
iei j=o is/ 

Thus by the induction hypothesis we have nA{n) > and therefore A(n) > 0. □ 

For our next result we will need the following lemma. 

Lemma 7. Exercise 22, p. 248] IfneN, then 

R(n) = 4er(n) - 4a(n/2) + 8a(n/4) - 32a(n/S) > 0. 
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Proof. If 8 | n, then 

R(n) = 4<j(n) - 4cr(n/2) + 8cr(n/4) > 4<r(n) - 4cr(n/2) > 0. 
If 8 | n, say n = 8k, then repeatedly application of identities ((T|) we obtain 
r(8Jfc) = 4cr(8fc) - 4a(4fc) + 8<r(2fc) - 32cr(fc) 
= 4o-°(8fc) + 4cr°(4fc) + 16a (2fc) 
> 0. 

This proves the lemma. □ 
Theorem 12. If 

oo oc 

[] (1 + x") 4 (l + x 2n f{l + x 4 ") 4 = a{n)x n , 



n=0 



then a(n) > for all n G No- 



Proof. By induction on n. Clearly a(0) = 1 > 0. Suppose now that the statement 
holds for j = 0, 1, . . . , n — 1. Note that 

oo oo 

II (1 + x") 4 (l + x 2n ) 2 (l + x in f = J] (1 - a;")- 4 (l - a; 2 ") 2 (l - a; 4 ")- 2 (l - x 8 ") 4 , 

n— 1 n— 1 

which by Theorem [TJ leads to 

na(n) = 4cr(n) - 2<j E (n) + 2<7 0A (n) - 4cr ,8(«) + 

n-l 

5^ a(j) (4a(n - j) - 2cr B (n - j) + 2a 0A (n - j) - 4a , 8 (n - j)) . 
i=i 

Further, by ([T|) and Lemma [7] we get 

4cr(n)-2cr B (n) + 2cTo,4(n)-4 < 7o,8(n) = 4<r(n) -4cr(n/2) + 8o-(n/4) -32cr(n/8) > 0. 
Then by the induction hypothesis we have na{n) > 0, and thus a(n) > 0. □ 

References 

[1] T. M. Apostol, Introduction to Analytic Number Theory, Undergraduate Texts in Mathemat- 
ics, Springer, 1 edition, 1976. 

[2] L. E. Dickson, History of the Theory of Numbers, Vols. I-III, Chelsea Publ. Co., New York, 
1952. 

[3] M. El Bachraoui, Inductive Formulas for some Arithmetic Functions, Int. Journal of Number 
Theory, To appear (available on arXiv. 1103.5227). 
(1995). 

[4] G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Oxford University 
Press, USA, 6th edition, 2008. 

[5] K. Ono, S. Robins, and P. T. Wahl, On representation of integers as sums of triangular 
numbers, Aequationes Mathematicae, Volume 50, (1995), 73-94. 

[6] N. Robbins, Some identities connecting partition functions to other number theoretic func- 
tions, Rocky Mountain Journal of Mathematics, Volume 29, Number 1, (1999), 335-345. 

[7] K. S. Williams, Number Theory in the Spirit of Liouville, Cambridge University Press, New 
York, First edition, 2011. 

Dept. Math. Sci, United Arab Emirates University, PO Box 17551, Al-Ain, UAE 
E-mail address: melbachraouiOuaeu. ac . ae 



